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We study the motion of a string in the background of the Schwarzschild black hole in 
AdS^ by applying the standard arsenal of dynamical systems. Our description of the phase 
space includes: the power spectrum, the largest Lyapunov exponent, Poincare sections 
C^ and basins of attractions. We find convincing evidence that the motion is chaotic. We 

discuss the implications of some of the quantities associated with chaotic systems for 
aspects of the gauge/gravity correspondence. In particular, we suggest some potential 
relevance for the information loss paradox. 



1 Introduction 

The AdS/CFT correspondence [H |2l |3], H] has been mainly used as a window into the physics of 
strongly coupled field theories but it also provides an approach to some gravity mysteries. For 
example, the question of black hole radiation can be, in principle, reformulated as a process in 
a unitary field theory. 

The fact that the quantum numbers of certain operators or states in field theory can be well 
described by the corresponding classical string is the idea at the heart of Regge trajectories 
where the hadronic relationship J ~ M^ is realized by a spinning string. This idea originated 
in the 1960's and certainly predates the AdS/CFT correspondence, it has found many concrete 
realizations in this new context. More generally, the AdS/CFT correspondence provides a 
dictionary that identifies states in string theory with operators in field theory. One of the most 
prominent examples is provided by the Berenstein-Maldacena-Nastase (BMN) operators. The 
BMN operators [5] can be described as a string moving at the speed of light in the large circle 
of S*^, the operator corresponding to the ground states is given by O'^ = [1 / yJN^)Ti Z^ . 
Another interesting class of operators which are nicely described as semiclassical strings in the 
AdS^ X S^ background are the Gubser-Klebanov-Polyakov (GKP) operators discussed in [6]. 
They are natural generalizations of twist-two operators in QCD and in the context of A/" = 4 
supersymmetric Yang-Mills they look like Tr$^V(a^ . . . V a^)^^ ■ A very important property of 
these operators is that their anomalous dimension can be computed using a simple classical 
calculation and yields a prediction for the result at strong coupling A — S* = {\f\/Ti) In S. 

These, and similar semiclassical states, have been the main window into the spectrum of 
string theory in AdS^ x S^ with Ramond-Ramond (RR) fiuxes and have given very valuable 
information about anomalous dimension of operators with large quantum numbers. This is 
particularly crucial in the absence of methods that lead to the full spectrum of string theory in 
backgrounds with RR fiuxes. 

One of the most important results of XIX century mechanics was the identification of the 
role of the full phase space for characterizing a system, rather than the role of single trajecto- 
ries. In view of the success of individual classical trajectories of strings in the gauge/gravity 
correspondence, it is only natural that we pose the question of the meaning of the full phase 
space. As soon as we pose this question we have to reckon with the fact that solving for the 
trajectories in phase space, that is, reducing them to quadratures is almost never possible. 
Integrable trajectories, in this sense, are a zero measure set. Even more, there are a number 
of reasons to expect the deterministic dynamics of strings in AdS backgrounds to be chaotic, 
where by chaotic we mean moving as if ruled by chance. Starting from the 1898 paper by 
Hadamard [7], geodesic fiows of point-like particles on manifolds of constant negative curva- 
ture have been found to satisfy the so-called Anosov (U) condition |S1IS]- Roughly speaking, a 
dynamical system satisfies the Anosov (U) condition if near an arbitrary fixed trajectory the 
behavior of the neighboring trajectories with respect to the fixed one is similar to the behavior 
of the trajectories close to a saddle. This particular kind of instability is typical of chaotic dy- 
namics. AdS is a space of constant negative curvature and it is not difficult to show that, as in 
the Hadamard example, the corresponding Jacobi equation P [ID] indicates that the deviation 
between neighboring geodesies of a point particle undergoes an exponential divergence. Even 



if AdS is not, rigorously speaking, compact, its causal structure is such that null and timehke 
geodesies going to infinity are refiected back, remaining effectively confined pOj • This way, the 
dynamics of point particles in AdS may be very complex. Furthermore, in a interesting paper 
[llj, Frolov and Larsen showed that, even for the Schwarzschild black hole in asymptotically 
flat space, if we consider the scattering of ring strings instead of point particles, the corre- 
sponding geodesic flow appears to be chaotic. Not surprisingly, in this paper we flnd that the 
problem of a ring string orbiting around the Schwarzschild black hole in asymptotically AdS^ 
space is highly sensitive to the initial conditions. This is a trademark signature of chaos. We 
support this flnding by using a standard set of indicators which includes the power spectrum, 
the Poincare sections, the largest Lyapunov exponent and the basins of attraction. 

Therefore, if chaos is likely to be a common feature of the dynamics of semiclassical strings 
in AdS backgrounds, it is relevant to wonder what it implies for the dual quantum fleld theory. 
Our aim in this paper is to take the flrst steps into expanding the dictionary of the AdS/CFT 
correspondence into the realm of chaotic systems. We propose some potential interpretation 
on the fleld theory side of the sensitivity to the initial conditions which could be typical of the 
chaotic motion of classical strings. As an application, we briefly discuss how this interpretation 
could be potentially relevant for the resolution of the Information Loss Paradox. 

The organization of the paper is as follows. In section |2] we discuss the set up including 
the embedding of the string in the Schwarzschild black hole spacetime. We also present the 
Hamiltonian description of the system. Section [3] is concerned with a discussion of the numerical 
methods we use and the standard analysis of the solution. Section |4] contains a study of various 
signatures that generically indicate chaotic behavior. In section [5] we present an argument 
for the possible identiflcation of the operators being described by the ring string conflguration 
discussed in this paper as generalizations of the GKP operators. We also make an attempt to 
interpreting a quantity typical of dynamical systems in terms of the dual fleld theory, namely, 
the largest Lyapunov exponent. We summarize our ideas in the last section. 

2 Ring String in Schwarzschild-AciS's 

Looking for an example of chaos in AdS, we will describe the dynamics of a particular kind of 
string moving in flve dimensions. Although the details will become clear later, in flgure [T] it is 
presented a simplifled picture of the embedding of the string in spacetime. We start with the 
Polyakov action 

C = --^y/^g'^'G^^daX^dtX", (2.1) 

zna 

where G^i, is the spacetime metric of the flxed background, X^ are the coordinates of the string, 

gab is the worldsheet metric, the indices a, b represent the coordinates on the worldsheet of the 

string which we denote as (r, a). The background metric G^^ is of the form: 

ds"" = -fde + ^+r^ (de^ + sin^ ed^^ + cos^ ed(j)^) , /(r) = 1 + ^ - ^, (2.2) 

corresponding to the Schwarzschild black hole in AdS^. In order to explore the conflguration 
sketched in flgure [l| we consider a circular string whose world sheet is parametrized by (r, a) 



Figure 1: Four dimensional sketch of the ring string dynamics. 
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and is described by the following embedding 

t = t(r), r = r{T), = ^(r), = 0(r), i; = aa, (2.3) 

where a is a useful constant that will track the role of winding or more fundamentally the poten- 
tial differences between strings and particles. In the conformal gauge, the Polyakov Lagrangian 
takes the form 

1 
'2'Ka' 

It is worth noting that the term proportional to a enters with the same sign as the temporal 
direction, which suggests adding some energy to the system with respect to the particle limit. 
The equations of motion are 

ffn , /' 
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e = — rO - sin 9 COS 9 U^ + a"^), r^ cos^ 9 (p = I, (2.5) 

where E and / are constants of motion naturally related to the energy and angular momentum 
of this configuration, dot represents derivative with respect to r and prime with respect to r. 
The nontrivial constraint from fixing the conformal gauge is 

G^, {doX'^doX" + diX^diX") = 0, (2.6) 

and in the above background takes the form 

- fi^ + y'^ + r'^ (O^ + cos^ 9 0^) + a^ r^ sin^ ^ = 0. (2.7) 

The first three terms in the above expression can be thought of as describing the motion of a 
point particle with proper time r moving on the space given by (t, r, 9, (p) with a varying action 
given by the last term. A more suggestive way to write the constraint is 

r' + r'f9' + -J-^^ + a' r'f sin^ 9 = E\ (2.8) 

r"^ cos^ 9 

In which case it looks like the conservation of energy for a classical mechanical system. 

Looking back at figure [l} we note that in this configuration, the radius of the string changes 

in time as 

i?(r)=r sin ^. (2.9) 

There are different modes of motion which can be characterized by the range of values of the 
radius r. The radius of the string, R{t) above, also contributes to the intuition of the different 
modes. For example, on one possible mode of motion, the string can approach the black hole 
and get scattered, crossing over it. Then, it can get pulled back by the action of the black hole 
gravity, and so on. In this mode both quantities r{T) and 9{t) oscillate in a way that R{t) is 
always greater than the Schwarzschild black hole. Along these lines, and recalling the results 
in [TT], we thus anticipate three modes of motion: i) the ring string remains forever crossing 
back and forward over the black hole ii) the ring string completes a number of these oscillations 
before collapsing into the black hole, ii) the ring string completes a number of these oscillations, 
before escaping to infinity. 



2.1 Hamiltonian Approach 



Let us consider the motion in the Hamiltonian formahsm, this is particularly appropriate for 
the treatment of various quantities in the context of dynamical systems. We start with the 



Lagrangian (2.4) 
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The canonical momenta are 
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The Hamiltonian is of the form 
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Let us consider the equations of motion following from our Hamiltonian. For the t equation we 
have 



Pt = 0, t = -na —. 



(2.13) 



For the radial coordinate r 
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The equation of motion for 6 is 
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(2.14) 
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Finally, for the (f) coordinate, we have 



P4> = 0, = vra' 



r^ cos^ 6' 



Note that the constraint (2.7) takes the simple form of 

if = 0. 



(2.15) 



(2.16) 



(2.17) 



The equations of motion for t and can be eliminated from the system in favor of the 
corresponding constants of integration. Summarizing, our reduced Hamiltonian system is 
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p0 = —na— — 5— ar sinfc'cosfc' (2.18) 
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subject to the constraint 

Hr = 0. (2.19) 

To develop our intuition further we can think of this constraint as we did from its Lagrangian 



analog (2.7), in which case we simply have a point particle with a varying action. In other 
words, we could think of the above system as describing a point particle whose world line does 
not allow to be affinely reparametrized. 

2.2 Conserved quantities 

Classical conserved quantities of motion correspond to quantum numbers of the corresponding 
operator or state in the field theory. The worldsheet current of spacetime energy-momentum 
carried by the string is 

P^ = -7^^G'^.5"X^ (2.20) 

Zna 

where a' is the string tension and the indices a and fi were explained at the beginning of the 
section. Note that, interestingly 

P\ = — , P"<p = —■ (2.21) 

2na'' ^ 2na' ^ ' 

These quantities are related to the conformal dimension and spin of the corresponding operator. 
The quantity that is particular to our configuration is 

P% = —r^ sin^ 6 = —R^, (2.22) 

^ 2na' 2na' ' ^ ' 

where R is the radius of the ring string. Thus, for this quantity we anticipate three possible 
asymptotic behaviors: i) oscillatory, ii) approaching zero or iii) approaching a constant finite 
value. 



3 Numerical treatment of the system 

Let us briefly discuss aspects of tfie numerics tliat we feel are relevant for the reader. To solve 



the system (2.18) we used a seventh-eight order continuous Runge-Kutta method which, thanks 
to its adaptive scheme, provides great control upon the output accuracy. This accuracy can 
be assessed at any time by evaluating constraint (2.19). As it is shown in figure pi for a finite 



interval of time At the numerical solver warrants that \Hr\ < 6, where 6 is the (absolute and 
relative) error tolerance of the solver. The smaller the value of 6, the larger the corresponding 
value of At. Since our aim is to detect chaos in this system, we will be typically considering the 
asymptotical behavior of the corresponding solutions which would imply the need of using small 



values of S. Note that using constraint (2.19) to dynamically calibrate the trajectory is not 



recommended in this setup because we expect this system to be sensitive to the initial conditions 
and there will always be an error in the calibration coming from numerical truncation, that 



is, from solving if,, = to some numerical order. Fortunately, as we will show in section |4.2 
(see fig Is] and the corresponding table), our numerical results for very large r are robust in the 
sense that they converge uniformly to given values as 6 is decreased. This warrants that the 
first figures of our asymptotical numerical results obtained with reasonable values of 6 coincide 
with those to be obtained with smaller 6. Our numerical experiments yield that to get six 
representative figures we need to fix 5 = 10~^^. 

These numerical experiments and the results presented in the figures in this paper were 



obtained setting, in system (2.18), Mw^ = 1, a = 1, 6 = 10, a' = l/vr, E = 10 and / = 10. 
These ad hoc values were chosen to yield a horizon radius near unity [rn = 0.995085) and, 
though a corresponding analysis must be done at some point, we expect our results here to be 
generic with regards to the values of the parameters. 

As our working example, in figure IS] we present a plot of a solution for t^t) obtained with 
the following set of initial conditions: {r(0) = 5,p^(0) = 2.670855,0(0) = 0,pe = 51.559927}. 
Though it looks regular, it is difficult to observe any definite pattern behind the amplitudes 
and frequencies of this oscillation. Moreover, in figure |4] we plotted the projection into the 
{r,pr) plane of the phase space of the phase curve corresponding to this trajectory. It seems to 
fill densely a given region in the phase space. This complex behavior is typical of deterministic 
chaotic systems. The aim of the following section is to show that this is, indeed, the case for a 
ring string in the Schwarzschild black hole in AdSr,. 

We end this section with some comments about our programs for the numerical computation 
of the dynamical indicators described below which were implemented in Maple following the 
recipes in references jT2l |T3l IH] . They were checked and tuned using the canonical examples of a 
system of two decoupled harmonic oscillators and the well known Henon-Heiles system, which is 
a nonlinear nonintegrable Hamiltonian system, whose parameters allow both, for quasi-periodic 
as well as for chaotic behavior. 

4 Signatures of chaos 

The naive definition of deterministic chaos is that of a process that appears to proceed accord- 
ing to chance even though its behavior is in fact determined by precise laws (no randomness 



Figure 2: The reduced Hamiltonian as a function of proper time r for different values of the error 
tolerance S of the numerical solver. 




Figure 3: A solution r{T) 
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Figure 4: Phase curve corresponding to the solution plotted in pi projected into the (r, p^) plane. 
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involved) [121 [131 E] . In the case of interest to us, we will see that even if we have system (2.18) 
which describes the geodesic motion of the ring string, the practical impossibility of exactly 
setting the initial conditions leads to the impossibility of predicting the final outcome for all 
the geodesies with initial conditions set inside a given ball in phase space, no matter how small 
the radius of this ball. 

A rigorous mathematical definition of chaos is still lacking. In this direction there have not 
been fundamental advances since the works by Brown and Chua [151 ttSj • In those papers they 
presented several examples leading to the conclusion that chaos is more a philosophical term 
than a mathematical one; no matter what definition is given for chaos, there is always some 
example of chaos which cannot be proven to be chaotic from the given definition. Therefore, 
here we focus on showing that our system has a number of properties that make it difficult to 
predict the large time behavior for any given geodesic when compared with the correspond- 
ing asymptotics of a neighboring geodesic at a given point of phase space. As was already 
mentioned, the interpretation of our results is done by comparison with the equivalent results 
obtained for a regular system of two decoupled harmonic oscillators and a quasi-periodic and 
a weakly chaotic configuration of the Henon-Heiles system. 

4.1 Invariant sets 

Since our emphasis is on the large time behavior for any given geodesic, it is important to 



determine the possible (positive) invariant sets of the phase space of system (2.18). For an 
invariant set with respect to the fiow, every phase curve starting on it remains on it for ever. 
Particularly important examples are fixed points, limit cycles and Kolmogorov-Arnol'd-Moser 
(KAM) tori. 

4.1.1 Fixed points 

First we look for fixed points. Let us, for simplicity, focus on the values of parameters we are 
using in this paper. Note that, for finite r, the right hand side of the equation for r in system 



(2.18) has real roots if and only if p^ = 0. Similarly, the right hand side of the equation for 9 
is zero if and only if its corresponding conjugate momentum is zero. In turn, setting to zero 
the right hand side of the equation for pg give us sin 6 times a positive term, which leads to 
6^: = nir, where n = 0, 1, 2, ... . Using all these values in the equation for p^ we obtain the 
algebraic equation 0.0099r^ — 0.02r^ -|- 0.02r^ + 2r^ — 1 = with only one real positive solution, 
r^, = 0.706009. So this fixed point is located behind the event horizon. The eigenvalues of 
the Jacobian evaluated at this point are (+80.598481i, -80.598481i, 14.190516, -14.190516), 
indicating that this equilibrium is non-hyperbolic, possibly of saddle type, but clearly is not an 
attractor, because to be an attractor all the eigenvalues should have negative real parts. 



Now, if we take the large r limit of the right hand sides of system ( 2. 18 ) , it can be determined 
that there is also a fixed point at {r^ = oo, Pr* = 0,9^ = mi,pg^:), for any finite pe*. If p^^ = 0, 
then for radial trajectories with pr > {pr < 0), this fixed point is an attractor (repellor). If 
Pe* 7^ 0, then at a large but finite value of r, the angle will eventually change, and as soon 
as this happen, it will result in a large angular momentum which can cause pr to change its 
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sign from positive to negative and the geodesic to start a bounce toward the center. It can 
be verified that in this fixed point the Hamiltonian constraint is also satisfied. The conclusion 
seems to be that the set of trajectories escaping to infinity has zero measure. 

4.1.2 Limit cycles and KAM tori 



Our reduced system (2.18) is Hamiltonian. If it were integrable, then it would have two first 
integrals in involution. Every level set of these two integrals is a 2-dimensional torus in the 
4-dimensional phase space. This torus is invariant. Phase curves are densely wrapped around 
it; generically both frequencies are functionally independent and they change from torus to 
torus. In general, not only both frequencies but also their ratio will change from torus to torus. 
Tori with only one frequency, i.e., limit cycles associated with asymptotically periodic motion, 
may also exist, but the probability of landing on such a torus by a random choice of an initial 
point in the phase space, is vanishingly small. The motion in the 4-dimensional phase space of 
an integrable system is generically quasi-periodic. The KAM theorem [T71 [9l [12] states that for 
nearly integrable systems (integrable systems plus sufficiently small conservative hamiltonian 
perturbations), most invariant tori are not destroyed, but are only slightly deformed, so that in 
the phase space of the perturbed system, too, there are invariant tori densely filled with phase 
curves winding around them quasi-periodically, with a number of independent frequencies equal 
to the number of degrees of freedom. These invariant tori form a majority in the sense that 
the measure of the complement of their union is small when the perturbation is small. If one 
continues to perturb a KAM torus, it reaches a stage where the nearby phase space becomes 
self-similar (has fractal structure). At this point the torus is critical, and any increase in the 
perturbation destroys it. There are still quasiperiodic orbits that exist beyond this point, but 
instead of tori they cover cantor sets called cantori. Thus, the transition to chaos in Hamilto- 
nian systems can be thought of as the destruction of invariant tori, and the creation of cantori. 
For weakly chaotic systems it is common to find tori and cantori coexisting in the phase space 

m 



Power spectrum 

In the absence of randomness, the power spectrum of a signal from an integrable system allows 
us to determine if it is periodic or quasi-periodic and which are the relevant frequencies. In the 
corresponding spectrogram these frequencies are visible as vertical bars with heights propor- 
tional to the square of the amplitude related to each frequency in a Fourier decomposition of 
the signal. For instance, for the two decoupled harmonic oscillators the corresponding diagram 
consists of two vertical lines. Adding randomness implies adding new spurious vertical bars and 
modifying the power of the original ones. When (white) noise overpowers the (quasi)periodic 
signal, we are left with a fiat spectrum over a defined frequency band. Signals from determinis- 
tic chaotic systems often are similar to those of noisy quasi-periodic systems, the power of the 
noise depending on how strong chaos is ^2i [iHl HI] ■ 

Figure [3] might give the false impression of periodicity. Indeed, it looks like two or three 
peaks appear at regular intervals of time r which would correspond to a dominant frequency. 
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However, in figure p\ we present the (symmetric) spectrogram for tlie signal r(T) plotted in 
figure |3| As can be seen, there are around twelve distinct peaks in frequencies. Most of these 
peaks have slightly blurred basis which suggests that the system is weakly chaotic. 

Poincare sections 

Since our phase space is four-dimensional, each energy level set is three-dimensional. We 
fix one such level set by setting Hr = 0. This three-dimensional manifold, fibered by two- 
dimensional tori, can be represented in ordinary three-dimensional space as a family of con- 
centric tori lying one inside another |9l [12]. Thus, if our system were integrable, for almost 
all initial conditions, a phase curve would densely fill an invariant 2-dimensional torus in the 
three-dimensional energy manifold. It is then easy to construct a two-dimensional plane in 
the three-dimensional energy level set transversally intersecting the two-dimensional tori of our 
family (in a family of concentric circles in the model in three-dimensional Euclidean space). 
This plane is a co-dimension 2 Poincare section. A phase curve beginning in such a plane 
returns to it after making a circuit around the torus. As a result we obtain a new point on the 
same circle (or any topologically equivalent curve) in which the torus intersects the plane. In 
this way there arises a mapping of the plane to itself, the Poincare map. With this tool, the 
transition to chaos in Hamiltonian systems can be seen as the destruction of the circle in the 
Poincare section, being substituted by a Cantor dust, the more diffuse, the stronger the chaos 
[121 HSl E]. In figure [6] we present a Poincare section at 6* = for the solution plotted in figure 
[3j As we can see, it confirms the result of the previous subsection of our system being weakly 
chaotic. Note also that we could take any point in this section as the initial conditions and we 
would obtain the same results presented in the figures of this paper. 

4.2 Largest Lyapunov exponent 

One of the trademark signatures of chaos is sensitive dependence on initial conditions, which 
means that for any point X E Ai, there is (at least) one point arbitrarily close to A4 that 
diverges from X, where Ai is an invariant subset of the phase space and the divergence needs 
not to be exponential [121 [ISl HI] ■ Let us consider two points in phase space, Xq and Xq-I- AXq, 



each of which will generate an orbit in that space according to the flow given by (2.18). These 
orbits are parametric functions of the proper time. If we use one of the orbits as reference orbit, 
then the separation between the two orbits will also be a function of time. In a neighborhood 
of stable fixed points or attracting limit cycles, this separation decreases asymptotically with 
time, while it usually increases in the neighborhood of unstable fixed points or of repelling limit 
cycles. For orbits winding around KAM tori this distance eventually settles down to a constant 
value. Thus, this separation is also a function of the location of the initial value and has the 
form AX{Xo,t). In cantori, the function AX(Xo,r) will behave erratically. It is thus useful 
to study the mean exponential rate of divergence of two initially close orbits using the formula 

A= hm hm -In^^^^. (4.1) 



13 



Figure 5: The power spectrum for the signal plotted in figJS 
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Figure 6: The Poincare section at 9 = for the signal plotted in figjS 
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This number, called the characteristic Lyapunov exponent, is the largest of a spectrum {Aj} 
which describes the evolution in time of a ball around an initial point in phase space. The 
Lyapunov spectrum was rigorously defined by Osedelec [H] in terms of the asymptotic behavior 
of the eigenvalues of the Jacobian of the vector field given by the right hand sides of the 
corresponding dynamical system. For Hamiltonian systems, where the volume of phase space 
is preserved by the flow, ^ Aj = 0. What is important for us is the question of whether A 
is positive as it will indicate that two trajectories that start out very close in phase space 
diverge in time. If the motion takes place in a bound domain, a positive (largest) Lyapunov 
exponent is a good indication that the system may be chaotic. Crucially for our study, it was 
already established that chaos, as characterized by a positive largest Lyapunov exponent, is 
coordinate invariant [12]. Moreover, recently this analysis was extended to include the class of 
transformations that do not preserve the boundedness of the orbits [20j. This may be relevant 
in the case of AdS. 

To calculate A we implemented the algorithm provided by Sprott [131 EI] • From the defini- 



tion (4.1) we see that A is an asymptotic quantity. To verify this numerically we should expect 



to observe that, as time r is increased, A settles down to oscillate around a given value. This 
is indeed observed in figure [7j A comment about the precision of this calculation is in order. 
As we have already mentioned, because of the cumulative numerical error, after some time the 
constrain Hr = is no longer satisfied within the error tolerance 6. The time interval where 
the constrain is satisfied is inversely proportional to 6. So, we have repeated the calculation 
with increased precision. The result is presented in the table and in figure [8] 



-In (Error Tolerance) 


Largest Lyapunov Exponent 
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0.0031875419244818369 


12 


0.0031833206867751884 


13 


0.0031823399315688655 


14 


0.0031823381278734457 


15 


0.0031823375459681767 


16 


0.0031823374248104873 


17 


0.0031823374063400799 


18 


0.0031823374086545868 


19 


0.0031823374085976368 


20 


0.0031823374085858968 


21 


0.0031823374085854369 


22 


0.0031823374085854015 


23 


0.0031823374085853993 


24 


0.0031823374085853991 



We notice a numerical convergence in the results. This trend gives us confidence that the 
value we computed is truly indicative of a positive largest Lyapunov exponent. The small but 
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Figure 7: The largest Lyapunov exponent as function of time for the signal plotted in figJ3 
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Figure 8: The largest Lyapunov exponent as function of 5 for the signal plotted in figlsl 
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positive value confirms that this is a case of weak chaos. 

4.3 Basins of attraction 

There is another signature of the unpredictabihty characteristic of chaotic systems which dis- 
plays the sensitivity to the initial conditions. This signature is clearly defined in systems with 
more than one possible time-asymptotic behavior. It is possible, and even common, that at fixed 
values of parameters of a nonlinear system, more than one kind of outcome may be obtained. 
For instance, as we have already anticipated, in our case we have three possible outcomes: the 
geodesic crosses the event horizon (the ring string collapses into the black hole), the geodesic 
radially escapes to infinity (the ring string goes to infinity while its radius reaches a limit finite 
value) and the geodesies settles down to an oscillatory motion (the ring string crosses back 
and forth over the black hole forever, with its radius getting stretched and shrinked chaoti- 
cally). The set of initial conditions (more precisely, the closure of this set) which eventually 
leads to each particular outcome is called its basin of attraction. We may plot these basins by 
assigning a different color to each outcome, and then coloring all initial conditions according 
to their outcome. If the boundaries between these outcomes are smooth, then the dynamics 
is regular. Conversely, if the boundaries are fractal, the dynamics is chaotic [22]. A fractal 
is a non-differentiable structure and so cannot be removed by any differentiable coordinate 
transformation. Thus, relevantly for us, a fractal basin boundary also provides an observer 
independent signature of chaos [23] . 

In figure |9] we plot the basins of attraction in a large subset of our phase space. All these 
initial conditions correspond to the values of 6 and r shown in the figure, Pr = and pe given 
by the constrain Hr = 0. Each figure contains 325 x 325 points. Red points correspond to 
the outcome of the geodesic crossing into the event horizon, blue points to escaping to infinity, 
while green points correspond to the oscillatory chaotic behavior as in figure [3j The condition 



for capture is simply r < r^. The conditions for escape were defined as in section 4.1.1 and 
correspond to r* > lOOrj^. For green points we consider trajectories that, after 5000 iterations, 
neither were captured nor escaped to infinity. 

There are some well defined red and green regions in figure [9j As we already noted, in our 
case, initial conditions yielding escape to infinity form a zero measure set. The red and green 
zones are divided by boundaries which look fuzzy, as if containing finer structure. To verify 
this, in figure [Tol we did a zoom of the region bounded by the blue rectangle in figure [9j We have 



found that, indeed, our description of figure [9] is also valid for the zoomed region. Furthermore, 



we did a zoom of the region bounded by the yellow rectangle in figure 10 , which is presented in 



figure 11 We see that finer and finer structure arises as we zoom into the boundaries. Moreover 



if, for instance, we zoom now the region bounded by the purple rectangle in figure 11 we will 



obtain something similar to figure [TOj This way we have verified the self-similar fine structure 
of the basins boundaries, a clear evidence for their fractality. Crucially for us, there are areas 
where an arbitrarily small change in initial conditions leads to a completely different outcome. 
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Figure 9: Basins of attraction. 



G.5 




20 



Figure 10: Basins of attraction. A zoom of the region bounded by the blue rectangle in figl9 
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Figure 11: Basins of attraction. A zoom of the region bounded by the yellow rectangle in fig 10 
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5 Field theory interpretation 

5.1 Identifying the corresponding operator 

In this section we argue that the operators being described by the ring string configuration 
discussed in this paper are generahzations of the Gubser-Klebanov-Polyakov (GKP) operators 
described in [B]. Let us, for a moment, assume that there is no black hole in our background, 
that is, that we are considering the ring string in pure AdS^ which correspond to M = in the 



metric described by equation (2.2). We will further introduce a new radial coordinate p of the 



form r = 6sinhp, where b is the radius of AdS^ as described in equation (2.2). The metric is 
then the standard AdS^ in global coordinates: 

ds"^ = - cosh^ pdt^ + dp^ + sinh^ pi^de"^ + cos^ Odcj)^ + sin^ Odilj'^). (5.1) 

We can consider the same embedding of the ring string in this background as described in 



equation (2.3). We end up with the following Nambu-Goto Lagrangian 



C = a sin6'sinhpY (t)^ cosh^p ~ P^ ~ sinh^p(^2 + cos^ 9{4>y). (5.2) 

The key conclusion we want to draw from this exercise is that the spacetime energy and angular 
momentum are 

, t sin 6' sinh p cosh^ p 

A ~ / dadr- 



S* ~ / dadr 



{iy COsh^ p - p2 _ (^2 + cos2 ^(0)2) ginh^ p 

sin 9 cos^ 9 sinh p sinh p 
(t)2cosh p — p^ — (^2 + cos^^((/))^) sinh p 



(5.3) 



These expressions are simple generalizations of those appearing in [6]. We thus, hope that 
our operators are generalizations of GKP operators. One key observation is in order. The 
derivatives in our cases are somehow exchanged, namely our coordinate p depends nontrivially 
on T while the dependence of p in GKP was nontrivial in a. We can speculate that the a ^ t 
switch is a sort of worldsheet duality which we loosely associate with S modular transformation 
on the string side. 

It was already clear in the case of GKP [H] that the precise operator dual to the string 
is hard to write down explicitly; in particular, its form depends on the parameters used: an- 
gular velocity or range of the coordinate p. More precisely, the string configuration in GKP 
interpolates between typical Regge states for small size of the spinning string and twist-two 
operators when the length of the string grows parametrically. We, therefore, will simply note 
that in some regime, the general form of the operators dual to the string configuration we are 
discussing should contain, as building blocks, operators of the form 

Oa,s ~ Tr$^(D^)^$^. (5.4) 
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More generally we can consider some insertion into this operator of the form 

Ot'' ~ Tr$^D(,,... /},,$-... D,,)<|.^ (5.5) 

The above operator, is not a singlet with respect to 5*0(6) but we are only concerned with the 
approximate general structure. One can even, attempt to implement momentum conservation 
along the string as done in |5j. At the expense of being overly speculative we conjecture that, in 
some regime of parameters, the ring string of our analysis is dual to an operator of the general 
form 

O r. Y.il-'^^t'- (5-6) 

n k=l 

In the context of the ring string in the Schwarzschild black hole background the situation is 
even more complicated. It is very likely that the state we are describing contains an averaging 
in the canonical ensemble over a large number of operators. We will not pursue this issue 
further in this paper. 

5.2 Largest Lyapunov Exponents and Poincare Recurrences 

A key aspect in the analysis of the largest Lyapunov exponent is the explicit assumption that 
one tracks the distance between two nearby trajectories at late times. Algorithmically, it means 
that the early time part of the trajectories is discarded form the computation. 

We find an interesting analogy with the fact that it has been long suggested that many as- 
pects of black hole physics require an understanding of late times evolution. The gauge/gravity 
correspondence naturally extends to the context of black hole physics. In principle, given this 
equivalence of a gravity theory (containing black holes) and a field theory, unitarity of field 
theories provides a possible resolution for black hole puzzles like information loss. Maldacena 
suggested to follow the detailed very long time properties of correlation functions in |24j . 

Let us recall the structure of correlation functions in generic field theories. Initial perturba- 
tion (insertion of operators) of a given thermal system will be damped by thermal dissipation as 
long as the time scale is too short to resolve possible gaps in the spectrum (Heisenberg times). 
For a unitary system with a discrete spectrum the perturbation is prevented from dying out. 
In fact, for any required precision there will exist a time for which the correlation functions 
returns to its initial value within the required precision [25] . 

Let us be a bit more precise, following the presentation, for example, of |26j. As mentioned 
before, Poincare recurrences are characteristic of finite bounded systems that evolve in a unitary 
way. For bounded system with discrete energy one considers an arbitrary operator A which 
evolves in time as A{t) = e**'^y4(0)e~**^. The corresponding matrix elements can be computed 
in the energy basis as 

Gs(t) = e-^(^) J2 |A,,fe*(^-^^)*. (5.7) 

Ei ^Ej ^E 

The object of study is the function GE{t) which is a correlation function. The key observation 
is that, if we define Heisenberg time as tn = l/ct', where ou = Ini^Ei — Ej), then this is the time 
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scale that reflects the discreteness of the spectrum. For t <^tH the spectrum is approximately 
continuous. 

If the matrix elements of some operator A in the energy basis have frequency with F, 
the correlator will decay with characteristic lifetime of order F~^: GE{t) Standard dissipative 
behavior in F~^ <^ t ^ tn- For t > tn most phases in GEif) would have completed a period 
and the function GEif) starts showing irregularities; it is a quasi-periodic function of time. 
Despite thermal damping, it returns arbitrarily close to the initial value over periods of the 
order of the recurrence time. 

In a very clear and pedagogical work Barbon and Rabinovici [26j clarified the importance 
of the distinction between discrete and continuous spectrum. In particular, they highlighted 
the fact that the presence of a horizon implies a continuous spectrum and thus, in that limit, 
Poincare recurrences should not be observed. 

What is the field theory object or operation dual to the computation of the Lyapunov 
exponent? We argue that the distance between two string worlsheets is related to the correlation 
function of the corresponding operators that the worldsheets are dual to. This intuition arises 
from the a similar situation with Wilson loops discussed in [27]. More generally, the correlator 
between two worldsheets is, in principle, a third worldsheet. In our case we believe that such 
worldsheet is approximated by the distance between the two original woldsheets. 

Now, by definition, the Lyapunov exponent describes the late time behavior of the distance 
between two trajectories. We therefore would like to translate this result as the late time 
behavior of a correlation between the two operators described by the worldsheets in question. 
Moreover, assuming that the late-time behavior of the correlators is dictated by the Poincare 
recurrences we conclude with the suggestion that the Poincare recurrence time tpR be related 
to the largest Lyapunov exponent A as 

tpR = \. (5.8) 

One last caveat is in order. The largest Lyapunov exponent is computed in phase space, 
therefore we essentially find a lower bound for the Poincare recurrence time. 

Late time correlation functions are the central object in the Information Loss paradox. The 
argument above ties Poincare recurrences to unitarity in the presence of a discrete spectrum. 
Thus, the relationship between chaos and its description in terms of the Lyapunov exponent in 
string theory might point to an underlying restoration of unitarity. 

As in the cases of |5] and [6], we believe the reason we are able to observed this, otherwise 
very small effect is due to the enhancing effect of large quantum numbers, in this case the 
spin S. Schematically we have: e~^ — )■ e~^ '^ ~ ^^(1); that is, what would have been a 
suppressed effect gets enhanced thanks to the presence of a large quantum number. 

It is worth mentioning the academic character of the above speculation. The problem with 
recurrence times is that typically they are so unimaginably long that they have no physical 
relevance. Think of the factors e"^^* as clocks, and let there be N of them. The question of 
recurrence is now the question of how long it takes for a certain configuration of clocks to 
reappear to within some angular accuracy A0. The recurrence time is {A(f)/27i)~'^ /u, where 
l/u is the average of I/ujj. Note that for just N = 10, 1/u = Isec, and A/2tt = 1/100 (an 
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accuracy of 1%), the recurrence time is lO^^sec. which is much longer than the age of the 
universe 1251. 



6 Conclusions 

In this paper we have presented a detailed analysis of the phase space structure of a ring string 
in the Schwarzschild black hole in AdS^ space. We have found convincing evidence in support 
of chaotic behavior. As is usually the case in bounded nonintegrable dynamical systems, we 
expect this dynamics to be generic. 

The general question that we are starting to address is: What is the interpretation, on 
the dual field theory side, of the main quantifiers of chaotic behavior? In other words, we try 
to extend the AdS/CFT dictionary to include entries related to chaotic dynamics. We have 
proposed that one of the key indicators of chaos - a positive largest Lyapunov exponent - 
be identified, on the dual field theoretic side, with the appropriate bound for the time scale of 
Poincare recurrences. Many other indicators of chaos remained to be discussed, for example, the 
interpretation of various fractal dimensions which are very characteristic of chaotic systems and 
for which we found evidence in the study of basins of attraction. Similarly, the interpretation 
of some of the concrete indicators we computed in this paper have not been addressed: the 
power spectrum, the Poincare sections and the basins of attraction. 

We finish by stating a number of open questions that seem particularly interesting to us, 
these questions are of a more technical nature and do not necessarily depend on the conceptual 
framework given by the AdS/CFT correspondence. The standard approach to dynamical sys- 
tems naturally specializes to the case of point particles. To our modest knowledge, not much is 
known or written about strings. There are substantial differences, most notably, the Virasoro 
constraints. It would be interesting to study and clarify the role of strings, or more generally 
extended objects (p-branes) using the arsenal of techniques available for dynamical systems. 

Integrability, understood as the ability to write the equations of motion in quadrature, is 
usually assumed to imply the absence of chaos in dynamical systems. This intuition, again, is 
rooted in the particle interpretation. In the case of motion of strings we are now talking about 
integrability of systems with infinitely many degrees of freedom. The string has infinitely many 
degrees of freedom as it can oscillate in infinitely many modes exp(mo") and integrability of 
the strings gives infinitely many conserved quantities. Is this as restrictive as in the case of a 
finite number of degrees of freedom where one can use the action/angle variables to completely 
characterize the system? There is a precedent to this question that might indicate that the 
answer is rather elaborate. Namely, the Toda system is integrable, however, it has been shown 
that any polynomial truncation of the Toda potential displays chaotic behavior. A behavior 
along these lines might occur in the case of the string. If we allow it to have all the excitations 
exp(mcr) turned on, it could behave very differently from situations in which we freeze some of 
the excitations by virtue of Ansatze similar to the one considered here. 

Another very interesting aspect is the role of the asymptotic conditions, that is, asymptot- 
ically flat versus asymptotically AdS [2H] • For example, in a previous study of a ring string in 
asymptotically Minkowski space time [TT] , the structure of the basins of attraction seem to be 
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different from the one obtained in our case. More generally, we would like to cleanly pinpoint 
the role of asymptotically AdS spaces. 
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